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ABSTRACT 

We employ contemporaneous optical, UV and X-ray observations from the XMM-Newton 
EPIC-pn and Optical Monitor (OM) archives to present, for the first time, simultaneous spec- 
tral energy distributions (SEDs) for the majority of the Peterson et al. (2004) reverberation 
mapped sample of active galactic nuclei (AGN). The raw data were reduced using the lat- 
est pipelines and are all analysed consistently. The virial mass estimates from Peterson et al. 
(2004) allow us to calculate Eddington ratios AEdd for the sample using the bolometric ac- 
cretion luminosities determined directly from the SEDs. We calculate hard X-ray bolometric 
corrections ^2-iokeV for the sample and confirm a trend for increasing bolometric correction 
with Eddington ratio proposed in previous studies. Our comparison with previous work on 
these objects suggests that the OM bandpass may be less susceptible to intrinsic reddening 
than the far-UV peak of the thermal disc spectrum in AGN, yielding larger bolometric cor- 
rections than previous work: ^2-iokeV ~ 15 — 30 for X^dd ^ 0.1, ^2-iokeV ~ 20 — 70 
for 0.1 < Asdd ^0-2 and ^2-iokeV ~ 70 — 150 for Asdd ^ 0.2, but part of this increase 
could be attributed to spectral complexity preventing accurate recovery of the intrinsic lu- 
minosity in some sources. Long-term optical-UV variability contributes a second-order, but 
significant change to the total bolometric luminosity when comparing multiple observations 
for individual objects. We also consider the effect of a recently proposed correction for radia- 
tion pressure when determining black hole masses with reverberation mapping, and find that 
the revised mass estimates do not significantly alter the range of bolometric corrections seen 
but may yield a narrower distribution of Eddington ratios. 
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1 INTRODUCTION 

It has long been a quest of multi-wavelength astronomy to better 
constrain the shape of the spectral energy distribution (SED) of ac- 
tive galactic nuclei (AGN), in order to better understand the nature 
of accretion onto supermassive black holes (SMBHs). The estab- 
lished paradigm for AGN includes a thermal accretion disc emit- 
ting in the optical and ultraviolet (UV); a corona above the accre- 
tion disc emitting in X-rays by inverse compton scattering of UV 
disc photons; and a dusty torus which absorbs part of the UV disc 
emission and re-emits it in the infrared (IR). The bulk of the accre- 
tion luminosity is visible in the optical-to-X-ray regime for most 
objects without heavy obscuration, and for such objects the repro- 
cessed IR should be excluded to avoid double-counting part of the 
emission. 

The seminal work of Elvis et al. (1994 ) has provided the 'tem- 
plate' AGN SED for many previous studies. The mean radio-loud 
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and radio-quiet SEDs presented in their study allow the calcula- 
tion of average parameters for AGN SEDs such as bolometric cor- 
rections, which are a key requirement when estimating bolometric 
luminosities or accretion rates for objects. However, the spread in 
parameters for individual AGN is large, and average SED parame- 
ters do not reflect this variation. |Marconi et al.| { |2004| and Hopkins 
|et al.| ( [2006| present new 'semi-empirical' SED templates which 
take into account the UV luminosity dependence of the optical-to- 
X-ray SED shape and propose a dependence of bolometric correc- 
tion on luminosity. One useful application of this luminosity depen- 
dent bolometric correction is to feed it into X-ray background syn- 
thesis models, which convert the observed energy density in X-rays 
from AGN into the supermassive black hole (SMBH) mass density. 
Such studies are extremely illuminating in that they allow us to also 
constrain the accretion efficiency of AGN in the process. However, 
the suggested luminosity dependence of bolometric correction does 
not seem to be borne out in the real AGN population. The work 
of Vasudevan & Fabian (2007 ) ('VF07' hereafter) presents obser- 
vationally determined AGN SEDs for 54 AGN for which optical, 



© 0000 RAS 



2 R.V Vasudevan and A. C. Fabian 



far-UV and X-ray data were avialable, and present the X-ray (2- 
lOkeV) bolometric corrections ^2-iokeV (£boi/£2-iokev) deter- 
mined from them. In particular, they use data from the Far Ultra- 
violet Spectroscopic Explorer (FUSE) which provides a window 
onto the peak of the accretion disc emission in the far-UV which 
has been historically difficult to constrain. Their work suggests no 
clear luminosity trend for bolometric corrections, and highlights 
many uncertainties in the determination of SEDs. Some such un- 
certainties include intrinsic reddening in the optical-UV part of the 
spectrum, variability in X-rays, lack of simultaneity in the data and 
uncertainties in the black hole (BH) mass estimates used to con- 
struct the SEDs. However, they also propose a dependence of bolo- 
metric correction on Eddington ratio (AEdd = ^boi/^Edd, where 
L Ed d = 1.38 x 10 38 (M BH /M o ) for a black h ole of mass Mbh), 
which may now be reinforced by some studies ( [Kelly et"aL| 2008 
IShemmer et al.|2008) . 

Most AGN SEDs presented in the literature have not been con- 
structed from contemporaneous broad-band data. Since the degree 
of X-ray variability is typically far greater than the optical and UV 
variability, this may not be a significant uncertainty; but certainly 
when observations in different wavebands are separated by many 
months or years, the accretion state of the AGN represented in the 
UV disc emission is not likely to be the same as that shown by 
the X-ray observations. Brocksopp et al. (2006 ) present a collec- 
tion of SEDs for 22 Palomar Green (PG) quasars with simultane- 
ous optical, UV and X-ray data from the XMM-Newton satellite's 
EPIC-pn X-ray and Optical Monitor (OM) instruments, with a view 
to eliminate uncertainties in long-term variability present in non- 
simultaneous data. A brief analysis of their SEDs in VF07 reveals 
that SED parameters determined from simultaneous SEDs may be 
significantly different from those from non- simultaneous ones; the 
comparison is made difficult since the bandpass of the OM falls 
just short of the canonical peak of the 'big blue bump' (BBB) of 
accretion disc emission at ~ 1000 A, which the FUSE observations 
do reach. However, using the simultaneous OM and EPIC-pn data 
is extremely valuable since it removes a large number of variabil- 
ity related uncertainties. The OM band-pass may also be less prone 
to the effects of intrinsic reddening in AGN, and may therefore be 
useful in extrapolating a relatively 'un-reddened' optical-UV SED 
shape for AGN. 

The emergence of a possible Eddington ratio dependence of 
bolometric correction may provide an interesting window onto the 
different modes of accretion at work in different AGN. The work of 
She mmer et al.| ( [2006] ) highlights a correlation between the hard X- 
ray photon index F and AEdd for AGN, reinforced in their later 
work (Shemmer et al. 2008). Additionally, the Eddington ratios 
seen in Galactic black holes (GBH) are often known to be corre- 
lated with particular types of accretion states: the low-flux, hard 
states generally exhibit low values of Asdd whereas high-flux, 
softer states are typically accompanied by higher values of AEdd 
(Remillard & McClintock 2006). The potential parallels between 
this and the Eddington ratio dependence of AGN SED shapes em- 
phasise the need for accurate determinations of the Eddington ra- 
tio, which require not only accurate SED data but also accurate 
estimates for the BH mass. There are many different methods for 
determining BH mass in AGN including direct measurements us- 
ing the host galaxy stellar kinematics, nuclear gas motions or the 
well known Mbh -bulge luminosity correlation (Marco ni & Hunt| 
[2003) . 

One well-established technique for estimating SMBH masses 
is reverberation mapping (RM), which has the advantage that it 
does not need the high angular resolution which the above meth- 



ods require. This is a proven and powerful technique which allows 
the general structure and kinematics of the broad line region (BLR) 
in AGN to be constrained. The foundations and assumptions be- 
hind the technique are discussed in detail in the review of Peterson 
|& Horne| ( |2004l ), but the essence of the approach involves moni- 
toring the response of optical broad line fluxes from the BLR to 
variations in the underlying (UV) continuum flux. The velocity dis- 
persion AV estimated from the broad line- widths is related to the 
size r of the BLR as r oc (Al/) -2 . We then require the system to 
be gravitationally bound, and since the gravity of the central black 
hole (BH) must be responsible for this, we arrive at the expression 
M BH = fr(AV) 2 /G for the black hole mass, where / is a scaling 
factor determined by the more detailed geometry and kinematics of 
the BLR. Values for / between ~ 3 and ~ 6 have been proposed 
previously, by requiring that the reverberation mass estimates obey 
the Mbh — cr relation (where a is the velocity dispersions of stars 
in the galaxy bulge), as done by Onken et al. ( 2004 ) and Marconi 
|et al.| ( [2008] ). The factor / is still somewhat uncertain however, and 
this propagates into uncertainties of a factor of ~ 3 in the virial 
mass estimates from RM. 

One of the most detailed reverberation mapping studies to date 
is that of Peterson et al. ('2004), who perform a thorough and con- 
sistent re-analysis of all previously available broad line reverbera- 
tion mapping data on a sample of 37 AGN, finding BH mass esti- 
mates for 35 of them (table 8 of their work). The mass estimates 
for two AGN in their sample, NGC 4593 and IC 4329A have ex- 
tremely large error bars, but in the case of NGC 4593 a more re- 
cent, better constrained RM mass estimate is available from jDen^] 
|ney et sl\] p006| ), calculated using a newer determination of the 
time lag between the BLR and continuum. SEDs for 21 of the AGN 
covered by Pet erson et al.| ( |2004| were presented in VF07 but using 
older mass estimates from diverse sources in the literature and non- 
contemporaneous data. In this work we refine the approach em- 
ployed in VF07 and present newly determined SEDs for a total of 
29 objects from the Peterson et al. ( 2004 ) sample. 

Since the RM technique requires that the optical-UV spec- 
tra are monitored for a sufficient time period to be able to detect 
changes in continuum emission, this necessitates that the sample 
is sufficiently bright in the optical and UV. We therefore miss out 
on optically faint classes of AGN such as obscured/type II AGN, 
but a detailed investigation of the accretion emission in obscured 
objects is beyond the scope of this study. We instead make the 
most of the |Peterson et al. (2004) sample (augmented with the RM 
mass for NGC 4593 from Denney et al. 2006 ) as providing a set 
of mass estimates uniformly determined by one method, and com- 
bine this with simultaneous optical, UV and X-ray data from the 
XMM-Newton archives as done by Brocksopp et al. (2006 ). This 
allows us to construct simultaneous SEDs and calculate Eddington 
fractions and bolometric corrections with significantly reduced sys- 
tematic uncertainties due to long-term variability and discrepancies 
between mass estimation techniques. This work represents the first 
such determination of simultaneous SEDs for AGN with consistent 
mass measurements. 



2 DATA SOURCES AND REDUCTION 

The simultaneous multi- wavelength data on the AGN in the |Peter-] 
| son et al.1 |2004 ) sample were obtained from the XMM pn and OM 
instruments (with the MOS data only being used in the initial stages 
for source detection purposes). Usable data were not available for 
five AGN: PG 0026+129, Mrk 817, P 1617+175, PG 1700+518 
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(for which no XMM data were available) and PG 0804+761 (for 
which no OM data was available in the one available XMM obser- 
vation). The original data files ('ODFs') for the observations were 
downloaded from the XMM-Newton Science Archive (XSAf] We 
generally follow the approach of Brocksopp et al. (2006) in reduc- 
ing the data. As part of our initial checks, we downloaded a few 
ODFs used in their work and ensured that the main results (such as 
photon indices and luminosities) were reproduced by our approach, 
within errors and allowing for changes in the XMM Science Anal- 
ysis System (SAS) pipelines. A list of the observation ID numbers 
for the observations used is given in Table [T] below, along with de- 
tails of the dates, observation times and window modes used for 
each observation. 

The X-ray data were reduced using the SAS v7.1.0 EPCHAIN 
and EMCHAIN for the pn and MOS instruments respectively. Auto- 
matic source-detection was performed using the EDETECT_CHAIN 
routines where possible; if this routine failed for a particular source, 
the source and background regions were identified manually. Cir- 
cular source regions (default radius 36 arc sec) were used to produce 
source event files and source spectra. The pn only was used for all 
subsequent parts of the analysis, and the pn event file was filtered as 
specified according to the standard guidelines (allowing only pho- 
tons with pattern 0-4 and quality flag 0). Source-free regions near 
the source - preferentially on the same CCD chip - were used to 
extract background event files and spectra. The datasets were in- 
spected for pile-up using the EPATPLOT routine and the 10-12 keV 
background light curves were inspected for flaring. Annular source 
regions were used in sources where pile-up was detected, and the 
TABGTIGEN tool was used to set a count-rate threshold from which 
to calculate the usable time intervals in each observation. The re- 
sults of these checks are also recorded in Table [T] Response matri- 
ces and auxiliary files were generated using the RMFGEN and ARF- 
GEN tools, and the final pn spectra were grouped with a minimum 
of 20 counts per bin using the GRPPHA tool. 

The Optical and UV data from the OM were reduced using 
the OMICHAIN pipeline. Point source and extended source iden- 
tification is automatically performed as part of this pipeline, and 
the point source corresponding to the nucleus was selected using 
the images generated for each waveband, ensuring that the identi- 
fied optical or UV counterpart was coincident with the source re- 
gion used for the X-ray spectrum. Since the reverberation-mapped 
sample will be relatively bright by necessity, we do not make any 
specific attempt to correct for host-galaxy contamination to the 
nuclear magnitudes obtained, and in most cases, the nucleus was 
clearly identifiable in the OM images. The optical-UV photometry 
results were extracted from the SWSRLI results files, and the mag- 
nitudes corrected for Galactic extinction using the extinction law of 
Cardelli et aT] fl989). Values for the Galactic extinction E(B - V) 
were obtained from the NASA/IPAC Infrared Science Archivqj 
The data points from the OM were converted to XSPEC spectrum 
files using the FLX2XSP utility, part of the FTOOLS package. 



3 CONSTRUCTION OF SEDS 

The XSPEC 1 2 package was used to fit a basic accretion disk and ab- 
sorbed power-law model (DISKPN+WABS(ZWABS(BKNPOWER))) 
to the reduced pn and OM data. The DISKPN model takes the pa- 
rameters T max (maximum temperature in the accretion disc), R in 



(inner radius of the accretion disc) and K (overall normalisation). 
The inner radius Ri n was frozen at 6.0 gravitational radii, and 
the temperature T max was left as a free parameter. The normalisa- 
tion for the DISKPN model is defined as K = M^ u cos(i) /Dl/3 4 , 
where Mbh is the black hole mass, i the inclination angle, Dl the 
luminosity distance and /3 the colour-to-effective temperature ra- 
tio. SMBH mass estimates were taken from Peters orTet al.|(2004) 
(or |Denney et al.||2006| in the case of NGC 4593), and luminos- 
ity distances from Ned Wright's Cosmology Calculator (assuming 
a cosmology of H = 71 kms _1 Mpc _1 , Q M = 0.27). We per- 
formed brief tests to ascertain the importance of inclination angle. 
Adopting typical values between 0° and 35° can produce changes 
in K of up to 18 per cent, but the total bolometric luminosity is 
a function of not only K but also T max and the spectral fit in the 
X-ray regime. Despite the 18 per cent change in K due to inclina- 
tion angle, this propagates to give changes of only ^6 per cent in 
the bolometric luminosities calculated from the SEDs, which was 
generally smaller than the statistical error in the luminosities them- 
selves. Therefore, a value of i = 0° was adopted in all cases, and j3 
was assumed to be Unity. The normalisation thus determined was 
kept frozen. 

The BKNPOWER model was employed to fit the X-ray data, 
with the low energy branch of the model set to have negligible 
values within the energy range of the BBB. This ensured that the 
power law was not extrapolated to give spurious fluxes at ener- 
gies below the thermal disc spectrum; this is a particular danger 
for sources with soft (r < 2) hard X-ray spectra (see Fig. [TJ. The 
break energy between the two regimes was set to be 4 times the 
DISKPN T max parameter. The WABS model was used to incorporate 
Galactic absorption in the X-ray regime, with the Galactic neutral 
Hydrogen column density, N^ al ^ frozen at the value provided by 
the NH ftool. The ZWABS component allowed for some element of 
absorption intrinsic to the source and the intrinsic column density 
was left as a free parameter. Bad channels in the X-ray spectrum 
were ignored, and only the 1-8 keV (observed frame) data were 
used for fitting the X-ray continuum, to avoid the soft excess and 
any features above 10 keV such as the reflection hump. Because of 
the simple approach adopted for modelling the continuum here, we 
do not use the values of Nh from the resultant fits for any further 
analysis, and suggest that for many objects with complex absorp- 
tion, the far more detailed fits available in the literature should be 
referred to for a more accurate determination of their intrinsic col- 
umn densities. 

We inspected the resultant fits, and in cases where this sim- 
ple model yielded a best fit photon index V smaller than 1.5, we 
fixed r at 1.5 in a revised version of the fit on the grounds that 
such low photon indices would be unphysical in canonical inverse- 
Compton scattering models for the coronal emission. This was 
done for the objects NGC 3227, NGC 3516, NGC 3783, NGC 
4151, PG 1411+442, Mrk 335 (1) and PG 1307+085. The choice 
of T = 1.5 is broadly consistent with the values found in more de- 
tailed studies on these particular observations (Gondoin et al. 2003 



for NGC 3227, |Turner et al.|2005| for NGC 3516 and |Reeves et al. 



2004 for NGC 3783). We present the resultant SEDs for 29 ob- 
jects with plausible fits to the data in Fig. [2] (IC 43 29 A did not 
yield a satisfactory fit in the optical-UV). We calculate unabsorbed 
luminosities in the 2-10keV, 0.00 1-1 OOkeV (bolometric) bands us- 
ing a modified version of the 'fluxerror.tcl' code provided on the 
HEASARC websit^Jfor this purpose. Monochromatic unabsorbed 



1 http://xmm.esac.esa.int/external/xmm_data_acc/xsa/index.shtml 

2 http://irsa.ipac.caltech.edu/applications/DUST/ 



http://heasarc.nasa.gov/xanadu/xspec/fluxerror.html 



© 0000 RAS, MNRAS 000, 000-000 



4 R. V. Vasudevan and A. C. Fabian 



Object 


Redshift 


XMM 


Observation 


Window mode 1 


Observation 


Usable percentage 


Source 


Pile-up 3 






Observation ID 


date 




time (ks) 


of observation time 


Variability 
(AI rms /(I» 

2 






n mini 
U.U33U1 


m <OQ/im m 
Ul jZo4U1U1 


onno no oa 
ZUU3-Uo-ZO 


o W 


1 00 Q 

133. o 


QC A 

yj.o 


n noco 
U.Uoj / 




ion i 


U.UjoIU 


no no to no m 
UZU3 /ZUZU1 


onn/i 1 n no 
ZUU4-1U-U© 


QWT 

b w 


on a 
/U.4 


QA O 

o4. / 


n nooo 
U.Uo33 




a t-v 1 on 


U.U3Z / 1 


m /iti am m 
U14 / 1VU1U1 


onno no o/i 
ZUU3-U5-Z4 


O W 


1 1 O 1 

1 1Z.1 


Q/l Q 


n nonn 
U.U/UU 




rairaii y 


U.U4 /UZ 


Ai ai n/inom 
U1U1U4UZU1 


onnn no n< 
ZUUU-U/-UJ 


rr 


oo n 
33. U 


QO A 


n oqoc 
U.3o3j 


\/ 

I 




U.UIoUj 


m ah aac\\ m 
U14 /44U1U1 


onno no aa 
ZUU3-U5-UO 


o W 


1 oa n 
130. U 


qa n 
oO.U 


n nooo 
U.U/Z / 




A/T„1, 1 1 A 

MrK 1 1 U 


n no con 
U.U3jZ9 


non 1 1 oa<ai 
UZU1 13UjU1 


onn/i 11 1 < 
ZUU4-1 1-1 j 


QWT 


/IO /I 

4 1 A 


O/I /I 

V4.4 


n nQA/i 
U.Uoo4 




MrK z /y (i) 


U.U3U4J 


nono/i onAAi 
U3UZ4oUoUl 


onn^ 11 in 
ZUUj-1 i-iy 


o w 


OQ O 

3o.Z 


qq n 
oo.U 


n 1 ooc 

u.iyzj 




MrK Z /y (Z) 


U.U3U4j 


nono/i OA/1A1 
U3UZ4oU4Ul 


onn< 11 1 < 
ZUUj-1 1-1 D 


5W 


CO o 

jy.o 


QO O 

o3. / 


n nooc 
U.Uo / j 




MrK Z/y (3) 


U.U3U4J 


nono/i oa<ai 
U3UZ4oUjU1 


onn< 11 1 n 
ZUUj-11-1 / 


o W 


CO Q 


QO C 
O / . J 


n 1 c/i q 
U.1j4o 




MrK 33 j (1) 


U.UZj /o 


aci aai nom 
UjIUUIU/UI 


onno no 1 n 
ZUU/-U/-1U 


LW 


OO 

ZZ.o 


91.7 


n oooo 
U.Z333 




MrK 33 J (Z) 


U.UZj /o 


A1 A1 c\ac\\ m 
U1U1U4U1U1 


onnn 1 o o< 
ZUUU-IZ-Zj 


rr 


OA o 
30. >> 


QQ O 
OO.Z 


n ooon 
U.3ZZU 


\/ 

I 


MrK j\Jy 


n no a ac\ 
U.U344U 


A'JAAAflA/IAI 

U3UOUVU4U1 


onnA c\a oc 
ZUUo-U4-Zj 


ib w 


on n 
/U.U 


OO o 

vz.3 


n noo/i 
U.U/94 




MrK j^U 


U.UZo3o 


nom nonom 
UZU1UZUZU1 


onn/i no c\a 
ZUU4-U/-U4 


o w 


1 1 o o 
1 1Z. / 


QA 1 

oO.l 


n conn 
U.jZUU 




MrK /y 


n mo 1 n 
U.UZZ19 


m AIO^AOAI 

UlU3ooUoUl 


onnn 1 n nn 
ZUUU-1U-U9 


ib w 


o c 

Z.J 


QQ A 
OO.O 


n oo 1 1 
U.ZZ1 1 




JNLrL. 3ZZ / (1; 


U.UU3oo 


A/1AA07A1 m 

U4UUZ /U1U1 


oaaa 1 o no 
ZUUO-1Z-U3 


LW 


1 no o 

iu/.y 


OO o 

yy.Z 


n 1 ooc 
U.1Z3J 




JNCjC 3ZZ / (Z) 


U.UU3oO 


m m a^aiai 
U1U1U4U3U1 


onnn 1 1 oo 
ZUUU-1 1-Zo 


FF 


/in 1 
4U.1 


OO 1 


n 1 oqo 
U.lVoZ 




MP P O C 1 A 

JNLrL. 3j1o 


U.UUoo4 


m A7/IAA7A1 

U1U /4oU /Ul 


onm 1 1 no 

zuui-i i-uy 


o W 


1 on n 
13U.U 


O/I o 

y4. y 


A 171 C 

U.l /1j 




MP P IICJ / 1 \ 

JNCjC 3 /o3 (1) 


n AAOT3 

U.UU9 /3 


m 1 OO 1 ACA1 

Ul IZZIUjUI 


onm 1 o 1 o 

zuui-iz-iy 


ib w 


1 oo o 
13 /.o 


OA /I 

9o.4 


n 1 onA 
U.13UO 




MPP llOO ZO\ 

JNLrL. 3 /o3 (Z) 


U.UUy /3 


m 1 oo 1 nom 
Ul 1ZZ1UZU1 


onm 10 n 
ZUU1-1Z-1 / 


QWT 


1 OO Q 

13 /.o 


OO Q 


n 1 /loc 
U.14ZJ 




MPP /IAC1 /1\ 

JNCjC 4UM (1) 


U.UUZ34 


m nm a 1 /im 
UlUy 1414U1 


onm ac 1 a 
ZUUI-Uj-Io 


ib W 


1 oo n 
1ZZ.U 


on a 
9U.0 


n OOAQ 

U.3900 




MPP /lfl<1 ZO\ 

JNLrL. 4Uj 1 \Z) 


U.UUZ34 


Ul j /joUlUl 


onno 1 1 oo 
ZUUZ-1 1-ZZ 


LW 


C 1 o 


OQ O 

9o.Z 


n oooc 
U. / /yj 


\/ 

I 


MP P/I1C1 

JNCjC 41M (1) 


U.UU33Z 


U143jUU1U1 


onno. a< oc 
ZUU3-Uj-Zj 


ib W 


1 o n 
19. U 


QA 1 
OO.l 


n nooo 
U.Uo39 




MPP A 1 < 1 ZO\ 

JNLrL. 41 J 1 (Z) 


U.UU33Z 


m /IKAAOA1 

U143jUUZU1 


onno n< oa 
ZUU3-Uj-Zo 


QWT 
b W 


ion 
lo.y 


OO o 

yy.Z 


n noo/i 
U.Uo /4 




MPP A 1 ^ T 

JNCjC 41 j 1 (3) 


U.UU33Z 


m /I O CAA1A1 

U143jUU3U1 


onno a? oo 
ZUU3-UJ-Z / 


ib W 


ion 

Lo.y 


OQ Q 
90.0 


n nooo 
U.U9Z9 




MPP /I ^O^* 


n aaoaa 

u.uuyuu 


uiuyy /U1U1 


onnn no no 
ZUUU-U/-UZ 


b W 


OQ 1 

Zo.l 


OQ A 

/0.4 


n 1 no a 
U.1U30 




MP P C C A Q 


U.U1 /I / 


UlUyyoUlUl 


onnn 1 o o/i 
ZUUU-1 Z-Z4 


ib W 


OA 1 

Zo.l 


OO 1 

vz.i 


n no ac 
U.UVoj 




MPP T/IAO 


U.U103Z 


m 1 o 1 om m 
Ul 1Z1 /U1U1 


onnn 1 o oa 
ZUUU-1 Z-Zo 


QWT 
b W 


1 o n 

iy.u 


O/I o 
94.3 


n 1 ono 
U.13U9 




dp nn^o iOci 
rLr UUjZ+Zj 1 


n 1 CCAA 

U.IjjUU 


nom a CA/1A1 
U3U14jU4U1 


onnc aa o/; 
ZUUj-Uo-Zo 


ib W 


on c 
ZU.J 


Q/l O 

o4.3 


n 1 ooc 
U.l /yj 




DP HO A A i 1AC\ 

ru Uo44+34y 


n aa /inn 
U.U04UU 


m noA^nom 
U1U300UZU1 


onnn 1 1 n/i 
ZUUU-1 1-U4 


rr 


OA A 

Z0.4 


QO A 
0^.0 


n /iooo 
U.4Z9 / 


\/ 

I 


DP AOC'J i A 1 A 

rLr U>03+414 


n oq /i 1 n 
U.Z341U 


ai i i onnom 

ui i izyuzui 


onm 11 oo 
ZUU1-1 1-ZZ 


LW 


1 C A 
lj.0 


OQ O 

9o. / 


n 1 CAQ 

U.Ijoo 




dp ion i i/ii / 1 \ 
rLr 1Z1 1 + 143 (1 ) 


n AOAAA 

U.UoU^U 


UZUoUZUlUl 


onn/i aa 01 
ZUU4-U0-Z1 


LW 


An n 
OU.U 


QO Q 

o3.o 


n ooco 
U.3Zj3 


\/ 

I 


PG 1211+143 (2) 


0.08090 


0112610101 


2001-06-15 


LW 


55.7 


84.2 


0.4841 


Y 


PG 1226+023 


0.15834 


0414190101 


2007-01-12 


sw 


78.6 


94.5 


0.0435 




PG 1229+204 


0.06301 


0301450201 


2005-07-09 


sw 


25.5 


100.0 


0.4125 




PG 1307+085 


0.15500 


0110950401 


2002-06-13 


LW 


14.0 


94.1 


0.2914 




PG 1411+442 


0.08960 


0103660101 


2002-07-10 


EF 


41.8 


85.1 


1.3188 




PG 1426+015 


0.08647 


0102040501 


2000-07-28 


FF 


17.6 


93.5 


0.3704 


Y 


PG 1613+658 


0.12900 


0102040601 


2001-04-13 


FF 


12.8 


67.5 


0.4279 




PG 2130+099 


0.06298 


0150470701 


2003-05-16 


SW 


37.9 


91.3 


0.2766 





Table 1. Details of XMM observations used in this paper. 1 - SW: Small Window mode, FF: Full Frame mode, LW: Large Window mode, EF: Extended Full 
Frame mode. 2 - Variability recorded in the rest frame 2-5 keV band. 3 - Sources marked 'Y' were corrected for pile-up as detailed in the text. * For NGC 
4593, the better constrained M-qh estimate from Denney et al. ( 2006 1 was used. Peterson et al. ( 2004) M-qh estimates were used for the rest of the sample. 



luminosities were determined at 2500A and 2 keV for determining 
the optical-to-X-ray spectral slope, aox- The calculated values for 
^2-iokeV, £boi, AEdd, ^2-iokeV and aox for each of these objects 
is presented in Table [2] 

4 DISCUSSION 

4.1 Checking the aox - (2500A) relation 

As in VF07, we check the SEDs generated against the 
aox — Lv(2500A) relation from the literature. It is known from 
the literature that the near-UV-to-X-ray spectral slope varies with 
the near-UV luminosity, and we expect our SEDs to be consis- 
tent with previous determinations of this relation. For definitions of 
aox and associated quantities, please refer to Strateva et al. ( 2005 ). 



The results are shown in Fig. [3] along with the recent best-fit deter- 
mination of this relation from Steffen et al. (2006 ). The calculated 
values of aox lie within the scatter expected from the |Steffen et al.| 
(2006) relation, implying that the SED shapes produced are in line 
with those expected from the literature. We note the presence of the 
X-ray weak quasar PG 1411+442, exhibiting a substantially lower 
aox than the rest of the population. 

4.2 Results for Individual AGN 

4.2.1 IC4329A 

Inspection of images from NED reveal that the host galaxy for this 
AGN is an edge-on spiral galaxy. As a result, the optical and UV 
continuum will be heavily contaminated with light from the galaxy 
and reddened by the pronounced dust lane visible in the optical 
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Figure 2. Spectral Energy Distributions for the AGN in the Peterson et al. ( 2004) sample. Filled circles represent the XMM pn and OM data, crosses represent 
the model fit with absorption and the dotted line represents the model with absorption removed. 



images. The OM points for this object do not produce a plausible 
fit to the DISKPN model at all, and since the mass estimate from 
Peter son et al.| ( |2004] ) has very large errors, we exclude this object 
from our SED sample altogether. 



4.2.2 NGC3227 

This AGN is classed as a Seyfert 1.5 in NED. As already men- 
tioned, it is also known to have significant optical and UV redden- 
ing ( |Crenshaw et al.||2001| > and our SED fit to the OM data con- 
firms this. For the observation designated NGC 3227 (2), the point 
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Figure [2] (continued) 



in the UVM2 band (2310A) is lower in vF v than the UVW1 point 
(29 10 A), contrary to the increasing shape expected from the UV 
disc emission, implying that the far UV has been significantly red- 
dened. The nucleus also seems to exhibit considerable variability 
in the UVW1 band, as the flux at this energy in observation (1) is a 
factor of ~ 2 larger than that for observation (2). VF07 adopt a de- 



redenned optical-UV continuum shape following Cren shaw et al.| 
(2001), but we leave the OM points as they are here. Additionally, 
this object is known to have significant and possibly variable X-ray 
absorption from Go ndoin et al.|(2003) , who analyse one of the same 
sets of XMM observations used here. Our simple ZWABS model for 
the intrinsic absorption yields values of log(A^/cm -2 ) between 
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Figure [2] (continued) 



^21.1-22.6, broadly consistent with the estimates of ^21.3-22.0 
from Go ndoin et al.| ([2003]). The presence of variable absorption 
presents difficulties in constraining the intrinsic X-ray luminosity. 
Our value of L2-iokeV ~ 3.2 x 10 41 erg s _1 is comparable to the 
value rsj 2.2 x 10 41 erg s _1 from the more detailed modelling of 
|Gondoin et al.| ( [2003] ). The two archival XMM observations used 



here were separated by six years, and the SED shape and derived 
parameters do not change appreciably between the two, as seen in 
Figs.|4]and[2] 
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Figure [2] (continued) 
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Table 2. Luminosities, accretion rates, hard X-ray bolometric corrections and cnox values for 29 AGN from the Peterson et al. ( 2004) sample, f - The revised 
mass estimate from Denney et al. (2006) was used for NGC 4593. Objects marked with an asterisk (*) had the photon index fixed at T = 1.5 (see ^3]for 
details). Errors are omitted for all quantities (except bolometric corrections) as they are generally small and systematic uncertainties are expected to dominate. 
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4.2.3 NGC3783 
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Figure 1. An illustration of the model combination used (wihout intrinsic 
or Galactic absorption components). The dotted line represents the diskpn 
model, the dashed line represents the bknpower model and the solid line 
shows the resultant total model. The bknpower model was employed to 
ensure that the X-ray power law was not extrapolated to UV and optical 
energies, below the Big Blue Bump energy range. The power law index for 
the low-energy branch of the broken power law was set to be - 1 . 
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Figure 3. X-ray loudness aox against (2500 A). Filled circles represent 
the a ox values determined from our sample; the black line is the best fit 
from[Steffen et al. ( 2006 1 and the grey shaded area shows the spread in their 
best fit. 



This object is known to have an outflowing warm abso rber { Kaspi 
|et al.||200l) , as discussed in detail by |Reeves et al.| p004| ), and 
we obtain a value of log(A^H/cm _2 )~ 21.3-21.4 from our sim- 
ple ZWABS model fit. This differs significantly from the value 21.7 
in |Reeves et al.| ( [2004] ), presumably because of their more detailed 
model fit to the data. These authors propose that the warm absorber 
is relatively stable between multiple observations, implying that 
variations in the SED shape are more likely to be due to intrin- 
sic continuum variation. The two observations used here (separated 
by two days) reveal very similar SED continuum shapes and the 
spread in bolometric corrections obtained is small (a factor ~ 1.2). 



4.2.4 NGC4151 

The host galaxy for this object is a barred spiral galaxy. It is also 
known to host high absorption and possesses a significant outflow 
( Krae mer et al.|2 005). This presents similar difficulties as discussed 
for the above objects for determining the intrinsic X-ray luminos- 
ity. We have processed three observations spanning three days for 
this object, which reveal a comparatively large spread in bolometric 
correction. Our simple ZWABS model is probably insufficient to ac- 
count for the complex nature of the outflowing absorption, and the 
underlying SED continuum could be more steady than suggested 
by the variable bolometric corrections. More detailed modelling of 
this AGN may yield a more constant SED shape. 



4.2.5 PG 1411+442 

This PG quasar is known to be X-ray weak as discussed by 
Brinkmann et al. (2004 ), with high intrinsic absorption (these au- 
thors find \og(Nn /cm"V 23.3). The bolometric correction de- 
termined for this object is ~ 1700, but is likely to be significantly 
overestimated if the intrinsic luminosity has not been recovered ac- 
curately. In VF07, another X-ray weak object was reported, PG 
1011-040, but this object did not show signs of heavy intrinsic 
absorption according to previous studies (Gallagher et al. 2001| ). 
Whether their low X-ray luminosities are intrinsic or not, the X- 
ray weak class of AGN, although small in number, continues to 
present a problem for the determination of a more universal AGN 
SED shape. 

4.2.6 NGC4051 

We present results for two observations for this Narrow Line 
Seyfert 1 (NLS1) AGN, separated by a year and a half. These two 
observations were analysed in detail by |Ponti et a l. (2006), who 
note the large difference in X-ray flux between the two observa- 
tions, distinctly harder spectrum for the low flux observation, and 
presence of a warm absorber in the spectrum. Other studies have 
sought to understand the source of the X-ray spectral variability 
seen in both high and low flux states (Uttley et al. 2004}, the de- 
tailed physics and geometry of the warm absorber seen in the long 
observation ( Krong old et al.|20 08 ]> and the complex emission and 
absorption line systems seen by |Ogle et al. (2004) in the long ob- 
servation. Mechanisms which have been suggested for the X-ray 
spectral variability in this source are variable absorption along the 
line of sight (partial covering), spectral pivoting, or the presence 
of a reflection component from the accretion disc. The analysis of 
Ponti et al. (2006) favours the last interpretation, and argues that 
the warm absorber, distinctly seen in the low-flux observation, is 
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present during both observations. If the large flux variations are due 
to underlying continuum variations, then the large spread we see in 
bolometric corrections (a factor ~ 5) could be typical for such ob- 
jects. However, it seems that recovering the intrinsic luminosity is 
sufficiently complex here that the SED shape is again difficult to 
constrain. 

4.2.7 Mrk 335 and PG 1211+143 

These objects are both NLS1 AGN. We present two observations 
for each of these objects, separated by almost seven years in the 
case of Mrk 335 and three years for PG 1211+143. The analysis 
of the observations used in this study and the previous analysis 
in VF07 both place them well into the 'high' Eddington regime 
(AEdd >^ 0.1), with X-ray spectra which resemble the 'soft' states 
of GBHs. The 20 ks observation of Mrk 335 is part of the Target 
of Opportunity programme in which XMM observations were ob- 
tained contemporaneously with Swift observations. This captures 
Mrk 335 at a remarkably low flux state ( |Grupe et aT[2 008 ). 

For high Eddington objects, it is likely that complex outflows 
may play a role in blurring our view of the nuclear emission. It 
is possible to model the spectrum with partial covering or strong 
reflection due to light bending ( Miniutti & Fabian 2004) for both of 
these objects, again making the intrinsic X-ray luminosity difficult 
to determine. This would naturally explain the large spread seen in 
the bolometric corrections (a factor ^2.9 for Mrk 335; a factor 
~ 3.7 for PG 1211+143). If the intrinsic nuclear X-ray luminosity 
is actually higher, this would tend to move the objects to higher 
Eddington ratios and smaller bolometric corrections. 

4.3 Revisiting bolometric corrections 

The simultaneous SEDs produced from our refined approach allow 
us to revisit the dependence of bolometric correction on Eddington 
ratio put forward in VF07, and the 2-10 keV bolometric corrections 
from the new SEDs are plotted against Eddington ratio in Fig. [4] 
(excluding the anomalous high bolometric correction for the X-ray 
weak AGN PG 1411+442). We note that nuclear emission from 
radio loud objects is likely to be contaminated with jet emission. 
We employ the values of radio loudness (R = Lq cu1 /Lb) presented 
in |Nelson| ( [2000t for this sample of AGN to identify three AGN (3C 
120, 3C 390.3, PG 1226+023 = 3C 273) with R > 10. 

We present the binned bolometric corrections against Edding- 
ton ratio with the radio loud objects removed in the left panel of Fig. 
[5] Error bars on the mean bolometric correction values for each bin 
are standard errors on the means obtained when averaging the bolo- 
metric corrections. For objects where multiple observations were 
available, the average bolometric correction and Eddington ratio 
calculated from all available observations were employed to cal- 
culate the bin average. It should be noted that in some cases, the 
dramatically different bolometric corrections seen for the same ob- 
ject in different observations could be reflective of an underlying 
change of state of the object, and in such a scenario a simple av- 
erage bolometric correction may not be meaningful. The effect of 
partial covering or light bending on the bolometric correction of an 
object are such that objects which intrinsically have values (A;, m) 
will appear at (Ao, fto) where Ao < \% and ko > k%. 

The right panel of Fig [5] shows the same results with the radio 
loud objects included, for comparison with Fig. 12 of VF07. With 
the consistent mass estimates and simultaneous SEDs presented 
here, we firstly confirm the previously identified trend for increas- 
ing bolometric correction with Eddington ratio. However there is 




Figure 6. Binned hard X-ray (2-10keV) bolometric correction against Ed- 
dington ratio for the sample with radio loud objects and low X-ray flux ob- 
servations (for NGC 4051, PG 1211+143 and Mrk 335) removed (red filled 
circles and lines). The VF07 results are presented for comparison (open 
grey squares and dotted grey lines). All other conventions are as specified 
forFig.[5] 



a marked shift towards larger bolometric corrections at all Edding- 
ton ratios. This could be because the band-pass of the XMM-OM 
is less significantly affected by intrinsic AGN reddening than the 
FUSE band-pass employed in VF07, giving larger bolometric lu- 
minosities. 

Our lowest Eddington bin could be subject to significant un- 
certainties, particularly because it includes NGC 3227 which is 
known to have significant reddening ( Crenshaw et al. 2001). How- 
ever, the increase in bolometric correction due to intrinsic redden- 
ing is likely to be as much as a factor of ~ 2 (see VF07) which 
would only increase the average for this bin by a few per cent from 
the uncorrected value of ^ 22. 

We note from Fig. [4] that for many of the objects with multi- 
ple observations where variable absorption or partial covering have 
been reported, there is a possibility that the low X-ray flux obser- 
vations are likely to give luminosities in the X-ray band signifi- 
cantly smaller than the intrinsic luminosity. The degree of variation 
in SED shape is particularly large for NGC 4051, PG 1211+143 
and Mrk 335 and so we plot bolometric corrections against Ed- 
dington ratio with the lower X-ray flux observations NGC 4051 
(1), PG 1211+143 (1) and Mrk 335 (1) excluded from the calcu- 
lations of average bolometric corrections (Fig. [6}. The very high 
average bolometric corrections seen in the high Eddinton ratio bins 
in Fig. [5] are now reduced substantially, reproducing the shape of 
the relation seen in VF07 more closely. 



4.4 Long term optical-UV variability 

The study of Smith & Vaughan ( 2007) demonstrates that the rms 
variability of the optical-UV continua in AGN is of the order of a 
few per cent, over the course of 40 - 350 ks observations. It is also 
instructive to compare the differences between optical-UV con- 
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Figure 4. Hard X-ray (2-10keV) bolometric correction against Eddington ratio for the sample. Black filled circles represent objects where only one observation 
was used. Empty circles represent the results from individual observations for objects with multiple observations. Blue filled squares represent the average 
values determined from multiple observations, and the blue 'bar & box' error bars show the ranges spanned by multiple measurements, or the intrinsic error in 
the measurement if this is greater. 



tinua over longer time periods. There exist extensive studies on con- 
tinuum and emission line variability in AGN by the 'International 
AGN watch' consortiurrj^] using International Ultraviolet Explorer 
(IUE), HST and ground based data. For example, observations of 
NGC 3783 presented by Reichert et al. ( 1994| show variations in 
flux of ~30 per cent at 2700A over an 8 month period, and even 
larger factors are seen over longer periods: a factor of ^2 in the 
5 340 A flux for Fairall 9 (observed over ^3 months - |Santos-Lleo 
et al.|l997t , and a factor of ^7 in 5100 the A continuum flux for 
NGC 5548 (observed over eight years - Peterson et al. 1999'). Here, 
we present the variations in OM photometry for the three objects 
with biggest spread in 2-10keV bolometric corection, namely Mrk 
335, NGC 4051 and PG 1211+143, to understand how this vari- 
ability might affect the optical-UV SED shape and luminosities 

4 http://www.astronomy.ohio-state.edu/ agnwatch/ 



calculated from the SED. We see from table [3] that changes of the 
order of ~ 20 per cent can introduce changes of up to a factor of 
~ 2 in the total UV disk flux, since small changes in flux in the 
optical-UV photometry will be extrapolated to larger changes over 
the whole energy range of the BBB. While the variation in hard X- 
ray bolometric correction for an individual object is still attributable 
to X-ray variability in the main, we see from this that the variation 
in the optical-UV continuum can be a very significant contributor 
as well. This supports the case for employing simultaneous X-ray 
and optical-UV data wherever possible in determining the total en- 
ergy budget of AGN. 



4.5 Comparison with VF07 

We note that 21 out of the 29 objects in this study were included in 
the previous study of VF07. This allows a comparison of the simul- 
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Figure 5. Left panel: Binned hard X-ray (2-10keV) bolometric correction against Eddington ratio for the sample; 5 objects per bin, 3 radio loud objects 
removed (see text for details). Right panel: binned bolometric corrections against Eddington ratio including radio loud objects (5 objects per bin, the highest 
Eddington ratio bin has only 3 objects), presented with VF07 results (open grey squares and dotted grey lines) for comparison. Error bars in bolometric 
correction are standard errors on the mean values of ft2-iokeV> calculated from averaging bolometric corrections in each bin. 
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Ara(V) 


Ara(B) 
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Am(UVWl) 
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disk ' disk 


Mrk 335 




0.090 


0.166 




0.248 




1.92 


NGC 4051 








0.227 


0.181 


0.200 


1.60 


PG 1211+143 






0.074 


0.0642 


0.038 


0.003 


1.03 



Table 3. Estimating the effect of long-term optical-UV variability on the total disk flux from diskpn (defined here as -Fo.ooi-o.ikev)- 



taneous SEDs with the non- simultaneous SEDs from VF07, where 
the UV data were taken in all but a few cases from FUSE obser- 
vations of Scott et al. (2004 ). Firstly, this allows us to re-assess the 
importance of simultaneous observations, and secondly, to under- 
stand how significant reddening may affect the SED shapes in the 
far-UV FUSE and XMM-OM bandpasses. We present a compari- 
son of a sample of SEDs from VF07 with the SEDs from this work 
in Fig. 17] and the hard X-ray bolometric corrections are compared 
in Fig. [8] We remove an erroneous optical SED point for the AGN 
NGC 3783, which was incorrectly named as NGC 3873 in VF07 
(this alteration does not alter the previously published bolometric 
correction). 

We note that the range of Eddington ratio and bolometric cor- 
rections obtained here is very similar to those obtained in VF07. 
At high Eddington ratio, the values for Eddington ratio from VF07 
tend to exceed those obtained in this study. This could be in part due 
to the masses being systematically smaller in VF07. The masses 
for the overlapping sources in VF07 were obtained from the previ- 
ous reverberation mapping catalogue of Kaspi et al. (2000) and the 
single-epoch virial mass measurements of Baskin & Laor (2005), 
along with a few objects with masses from Peterson et al. I p004] > 
(NGC 3227, NGC 3516, NGC 4593, Mrk 279). Eddington ratios 
for individual objects in this study range from a factors of ~ 0.1 to 
~ 2.6 of their values in VF07, and on average are a factor ^0.8 
of the VF07 values. The Pet erson et al.| ( |2004| ) masses are gen- 



erally larger than the values used in VF07, but the differences in 
Eddington ratio between the two studies can be ascribed equally to 
differences in black hole mass estimates and bolometric luminosi- 
ties. Bolometric corrections from this study range from factors of 
~ 0.5 to ~ 4.0 of the VF07 values, and on average are slightly 
larger in this study (by a factor ~ 1.34). We note that if NGC 3227 
is excluded from this average (because the obvious presence of red- 
dening in the spectrum used here prevents a meaningful compari- 
son with the reddening-corrected spectrum used in VF07), we find 
that the bolometric corrections in this work are on average around 
~ 1.38 times larger than the VF07 values. 



4.6 Average SEDs for high and low Eddington ratios, 
comparisons with a real accretion disc SED 

We present the rest-frame average SEDs (normalised at leV) for 
the AGN in the highest and lowest Eddington ratio bins in order 
to compare the SED shape characteristics (Fig. [9}. We adopt the 
higher X-ray flux observations where multiple observations were 
available, for the reasons outlined in section [43] The highest Ed- 
dington ratio bin has (AEdd) = 0.61 and we note the presence of a 
substantially larger accretion disc component, coupled with a softer 
X-ray spectrum than the average SED for the low Eddington ratio 
bin ((AEdd) = 0.01). This confirms the trends reported in VF07 
for the variation in SED shape with Eddington ratio. The average 
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Figure 7. Comparisons of SEDs from this work (data: filled black circles, model: black crosses) with a selection from VF07 (data: open red squares, model: 
red dots). 



high and low Eddington ratio SEDs have markedly different frac- 
tions of ionizing UV luminosity (Xi3.6-iooev/£boi)- The average 
low Eddington ratio SED in this study has Li3.6-iooev/Lboi ~ 21 
per cent, rising to ~ 59 per cent for high Eddington ratio (taking 
Lboi = ^o.ooi-iookev)- The average high and low Eddington ratio 
SEDs in VF07 display more pronounced differences; at low Ed- 
dington ratios the fraction is typically ~ 9 per cent, rising to ~ 54 
per cent for high Eddington ratios; but the less dramatic change 
seen in this work is expected because the range of Eddington ratios 
spanned here is smaller than that in VF07 (the central Eddington 
ratios of the high and low Eddington bins in VF07 span a factor of 
~ 140; here they span a factor of only ~ 50). 

In this work, we adopt a simple multicolour disc model 
(DISKPN) for the BBB. It is instructive to compare how the optical- 
UV data compare with more sophisticated accretion disc mod- 
els. |Ross et aT] ( |1992| ) present self-consistent calculations for 
Hydrogen-Helium accretion disc spectra in AGN, based on the 
standard a-disc model ( Shakura & Syunyaev 1973), including fea- 
tures such as the H-I Lyman edge and He-II Lyman a and j3 
emission lines. We reproduce the disc SED for an accretion disc 
surrounding a Mbh = 10 8 Mq black hole with Eddington ratio 
AEdd = 0.167, along with the equivalent black-body disc spec- 
trum for comparison. (Fig. Qo| reproduced from Figure 4b of their 
paper, with the authors' consent). The 'real' disc spectrum follows 
the equivalent black-body spectrum at low energies (E < 4 eV) 
but deviates at higher energies, where the H-I and He-II features 



are evident. The descending high-energy tail is shifted significantly 
from that for the equivalent black-body disc spectrum. However, 
since the H-I and He-II features are in the unobservable extreme- 
UV, simple black-body model used here is adequate for our pur- 
poses, since only near-UV data are available to us. The difference 
in total flux between the 'real' disc model and the blackbody is neg- 
ligible (of the order of ~ 5 per cent) and would not significantly al- 
ter bolometric luminosities. For a handful of sources in this work, 
the UV points are close to the turnover in a standard black-body 
disc spectrum. In these cases, the data points would lie at most a 
factor of 2 above the equivalent 'real' disc spectrum, since the 
model values differ by approximately this factor at the turnover. If a 
'real' disc spectrum were fit to these points successfully, the result- 
ing disc luminosities would increase by a comparable amount. This 
could only be significant for perhaps 6 objects (3C 390.3, Fairall 
9, NGC 3516, PG 1307+085, PG 1426+015, PG 2130+099) out of 
the 29 in the sample. We exclude NGC 3227 from this list because 
the low UV flux is already attributed to reddening. 



5 SYSTEMATICS 

5.1 Biases in the sample 

We note the possibility that the increase in bolometric correction 
with Eddington ratio could be influenced by the presence of Lboi 
on both axes. We assess the importance of this by generating a ran- 
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Figure 8. Comparison of bolometric corrections against Eddington ratios 
from the current work with the results from VF07. Black filled circles rep- 
resent results from this work; blue empty squares represent results from 
VF07. Error bars are omitted for clarity. For the typical size of error bars, 
refer to Fig. [4] 




E/keV 

Figure 9. Average SEDs for AGN in the highest ((A Edd > = 0.61, thick 
blue dashed line) and lowest ((Asdd) — 0.01, thick black solid line) Ed- 
dington ratio bins. The individual SEDs for AGN in those bins are depicted 
using thin blue dashed lines (( AEdd) =0.61 bin), and thin black solid lines 
«A E dd) = 0.01 bin). 




E / kcV 



Figure 10. Comparison between a more realistic accretion disc SED with 
the equivalent black-body spectrum. The curves are taken from Fig. 4(b) 
from Ross et al. (1992) with the authors' permission. 



dom sample of AGN with X-ray luminosities, redshifts and masses 
randomly chosen within the ranges shown in this sample. A conser- 
vative flux limit of 10 -15 ergs -1 cm -2 was employed in choosing 
luminosities. A random bolometric correction is assigned to each 
AGN between 5 and 1000, ensuring that the resulting bolometric 
luminosity does not lie outside the limits in our sample. A plot of 
bolometric correction against Eddington ratio calculated from these 
assumptions will be constrained to lie within certain limits, in par- 
ticular restricting the allowed bolometric corrections to lower val- 
ues at very low Eddington ratios. However, the range of scatter seen 
in the points is far greater than that seen in the data, implying that 
the limits imposed by plotting L ho i/L 2 -iokeV against Lboi/^Edd 
are not relevant here. The gradient implied from the random sam- 
ple is also steeper than that seen in the data. If a correlation is in- 
troduced between Mbh and L 2 -iokeV, me limits imposed by the 
data are narrowed somewhat, but such a correlation has not been 
reported, either in local or deeper samples (Pellegrini 2005 Babic 
|et al.|2007] >. These considerations suggest that the changes in SED 
shape/bolometric correction with Eddington ratio are intrinsic. 



5.2 Corrections to reverberation mapping masses due to 
radiation pressure 

Reverberation mapping masses have previously known to have an 
uncertainty of a factor ~ 3. However if this is symmetric, it is un- 
likely to change the trend of the bolometric correction-Eddington 
ratio plot significantly. However, recent work by Mar coni et al.| 
(2008) discusses the need to take into account the effect of radi- 
ation pressure when calculating virial mass estimates. They argue 
that in calculating the size of the broad line region (BLR) required 
for virial BH mass estimates, the effect of radiation pressure on the 
BLR clouds must be included as well as gravity. This becomes es- 
pecially important in AGN where the assumption that the system 
is gravitationally bound is likely to be violated, such as Narrow 
Line Seyfert 1 AGN (NLSls) which are thought to be accreting at 



© 0000 RAS, MNRAS 000, 000-000 



Simultaneous Snapshots of AGN from XMM 15 



high, possibly super-Eddington ratios. Their refined expression for 
the BH mass includes the standard virial term with the geometry- 
dependent scale factor /, and an additional term dependent on the 
radiation pressure scaled by a new factor g, which is dependent 
on the density of the BLR clouds and the SED shape of the illu- 
minating central source. They provide the empirically calibrated 
correction for radiation pressure in equation (6) of their paper, as a 
function of the 5100A monochromatic luminosity (ALa(5100A)). 

We try to assess the importance of this effect on Eddington 
ratios and bolometric corrections by applying a first order cor- 
rection to the reverberation mass values using the 5100 A lumi- 
nosities from our AGN SEDs. We adopt the values / = 3.1 and 
log(g) = 7.6 as found by Marconi et al. (12008 ) from requiring 
that the reverberation masses lie on the Mbh — cr relation. The 
average value of MggVMgg is - 1.35 for the whole sample. 



Marconi et al.| (2008> calculate ALa(5100A) more accurately by 
subtracting the host galaxy following the approach of |Bentz et al.| 
(2006), and find a similar average value for the whole sample, 
M B ^ /Mbh ~ 1.25 (obtained via private communication with the 
author). We do not find a systematically higher value for the four 
known NLSls in the sample (NGC 4051, Mrk 335, PG 1211+143 
and PG 2130+099), and some of the low Eddington ratio AGN have 
smaller masses with the new calculation, because of the lower value 
of / used ( [Peterson et al.|2004| use / = 5.5). 

We plot the bolometric corrections against Eddington ratio for 
the objects with individually corrected masses in Fig.[TT] eliminat- 
ing one object where the new mass value (and resultant DISKPN 
normalisation) results in an unphysical fit. The Eddington rates and 
bolometric corrections before correction are plotted for compari- 
son, and as expected, since the high Eddington ratio objects are 
those in which radiation pressure is likely to be higher, their cor- 
rected masses result in Eddington ratios lowered by up to ~ 1 dex 
with bolometric corrections changing by ~ 0.1 — 0.4 dex. Some of 
the revised mass estimates for low Eddington ratio AGN are now 
smaller, causing the lower Eddington ratio objects to shift towards 
higher ratios, revealing what may be an intrinsically narrower dis- 
tribution of Eddington ratios in the sample. Overall, the range in 
bolometric corrections seen is not altered significantly. 




Figure 11. Hard X-ray (2-10keV) bolometric corrections against Edding- 
ton ratio with masses corrected for radiation pressure according to |MarconT| 
|et al.| {2008) (filled circles), compared to results using the original reverber- 
ation mapping Mbh estimates (empty circles). 



on average for these AGN, which all seem to be in a low-Eddington 
fraction regime (0.008 < \Edd < 0.170). This would support 
the hints from VF07 and this work that ^2-iokeV is low for low 
Eddington fraction AGN. The distinctive accretion and SED char- 
acteristics of low-Eddington fraction AGN could be clarified by a 
similar analysis of a much larger sample of AGN using both X- 
ray and IR data. We would then be provided with a window onto 
the accretion modes at work in both 'unobscured' and 'obscured' 
AGN. 



5.3 Probing the low Eddington ratio regime 

It is clear that this sample contains few sources with Xsdd <0.1, 
apart from NGC 3227. This was also identified as a problem with 
the FUSE sample used in VF07, where supplementary UV data 
were gathered from other sources in the literature to augment the 
low luminosity/Eddington ratio regime. The tendency for redden- 
ing to be significant in this regime also presents problems. Addi- 
tionally, it is known from Fabian et al. (2008) that the majority of 
sources with significant absorption tend to cluster at low Eddington 
fractions. These observations suggest that it is necessary to use the 
re-processed IR to understand the accretion emission in low lumi- 
nosity Seyferts and AGN. Pozzi et al. ( 2007) present precisely this 
type of analysis for eight high-redshift luminous obscured quasars, 
combining XMM-Newton observations with Spitzer IRAC, MIPS 
24/xra and i^s-band observations. They remark on the likelihood 
that a significant proportion of luminous obscured AGN may es- 
cape identification due to faint optical counterparts, implying large 
obscuration and directing us naturally to the IR to study the re- 
processed nuclear (torus) emission. They estimate the bolometric 
luminosity by adding the total X-ray luminosity (0.5-500 keV) to 
the total IR luminosity (1-1000/im) and thus calculate bolometric 
corrections for their sample. Interestingly, they find ^2-iokeV ~ 25 



6 SUMMARY AND CONCLUSIONS 

We have presented, for the first time, a set of simultaneous optical- 
to-X-ray SEDs for 29 AGN from the |Peterson et"aT] {2004} 
reverberation-mapped sample. By using simultaneous data and 
consistent mass estimates, this work represents a refinement over 
VF07, yet reinforces their main conclusions. 

• In particular, we confirm an increase in bolometric correction 
with Eddington ratio, with bolometric corrections of ft2-iokeV ~ 
15 - 30 for A Edd < 0.1, ^ 2 -iokeV ^20-70 for 0.1 < A Edd < 
0.2 and ^2-iokeV ~ 70 — 150 for AEdd > 0.2, excluding three 
significantly (R > 10) radio loud objects. We note that the bolo- 
metric corrections from this study are generally higher than those 
in VF07 by a factor of ~ 1.3 — 1.4, and Eddington ratios tend to 
be marginally smaller in this study. The change in bolometric cor- 
rection is likely due to the reduced intrinsic reddening in the OM 
bandpass; this has the effect of increasing the overall optical-UV 
luminosity, giving higher bolometric corrections. Here we demon- 
strate that the changes in the SED with accretion rate AEdd are con- 
sistent with their results when performed on a sample with fully 
re-reduced simultaneous data and consistent mass estimates. 

• The availability of multiple observations for individual objects 
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pinpoints cases where constraining the intrinsic SED shape is prob- 
lematic because of the complex and variable nature of intrinsic ab- 
sorption, or the presence of outflows. NGC 3227 shows significant 
UV reddening in one of the available observations, but little over- 
all variation in the SED shape aside from this. NGC 3783 and NGC 
4151 both exhibit outflows, and for the latter in particular, this man- 
ifests as a relatively large srpead in bolometric correction (a factor 
of ~ 2). PG 141 1+442 is X-ray weak, but the large bolometric cor- 
rection of ~ 1700 is likely to be severely overestimated because 
of significant absorption. Conversely, it is also known that the X- 
ray spectra of X-ray weak AGN are not all accountable for with 
large amounts of absorption, and so they remain a somewhat un- 
known quantity in the larger 'puzzle' of the AGN SED. The two 
observations of NGC 405 1 in the XMM archive are known to have 
radically different X-ray fluxes and yield very different bolometric 
corrections, spanning a factor greater than ~ 5. A warm absorber, 
partial covering or a reflection component from the accretion disc 
could account for this large variation. Mrk 335 and PG 1211+143 
are both NLSls with high Eddington fractions, where the presence 
of outflows may confuse our view of the accretion emission. These 
considerations highlight that recovering the underlying continuum 
in AGN is a messy process and often requires a careful, object- 
specific approach. 

• Relatively small changes in the optical-UV photometry over 
months to years can lead to comparatively more significant changes 
in the total extrapolated flux in the BBB energy range, reinforc- 
ing the preference for simultaneous X-ray and optical-UV obser- 
vations where possible when calculating the total accretion energy 
budget. 

• Reverberation mapping provides a powerful and reliable 
method for getting SMBH mass estimates, which are a prerequisite 
for accurate Eddington ratio estimates. However, the uncertainty in 
the precise geometry and kinematics of the BLR contributes a fac- 
tor of no 3 uncertainty to the resultant mass estimates, but recent 
work by Marc oni et al.| (|2008) has highlighted a possible system- 
atic shift towards larger RM masses if radiation pressure is taken 
into account. We estimate that the ratio of the new to old masses 
for the Peterso rTet al.| ([2004) sample is ~ 1.35 on average, but 
with some low Eddington AGN actually exhibiting a decrease in 
mass. This does not significanlty alter the spread of bolometric cor- 
rections but may narrow the distribution of Eddington ratios in the 
sample overall. 

• It is clear that the low Eddington ratio regime is difficult to 
constrain with relatively bright Seyfert 1 AGN and quasars. The 
study of Fabian et al. (2008) on deep samples such as the Chan- 
dra Deep Field-South and Lockman hole suggests that the major- 
ity of sources are of low Eddington ratios. Moreover, they are often 
significantly obscured, and therefore optical-UV monitoring of the 
accretion disc emission is likely to be plagued by significant red- 
dening problems. The work of Pozzi et al. (2007) suggests that the 
bolometric corrections for low Eddington ratio AGN are likely to 
be small as obtained in VF07 and in this work, and a large sample 
of AGN with X-ray and IR data would then be useful for complet- 
ing the picture of AGN accretion emission across a much larger 
range of Eddington ratios. 

• We have demonstrated that good quality simultaneous optical, 
UV and X-ray data, coupled with consistent mass measurements 
continue to provide opportunities for understanding the accretion 
characteristics of AGN in new ways, particularly for the brighter, 
unobscured objects. The large archives from XMM-Newton and 
Swift provide great potential for further illuminating study. If such 
work is coupled with a consistent, well calibrated mass estimation 



technique, we will also be able to calculate accretion rates for the 
AGN from accurate determinations of their bolometric luminos- 
ity. Single-epoch virial mass measurements and various forms of 
the Mbh — Lhuige relation offer two such avenues, and a large 
database of SEDs with accretion rates for a well-defined sample of 
AGN remains an interesting subject for future work. 
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